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ABSTRACT 

The light curve of an exoplanetary transit can be used to estimate the plan- 
etary radius and other parameters of interest. Because accurate parameter esti- 
mation is a non-analytic and computationally intensive problem, it is often useful 
to have analytic approximations for the parameters as well as their uncertainties 
and covariances. Here we give such formulas, for the case of an exoplanet transit- 
ing a star with a uniform brightness distribution. We also assess the advantages 
of some relatively uncorrelated parameter sets for fitting actual data. When limb 
darkening is significant, our parameter sets are still useful, although our analytic 
formulas underpredict the covariances and uncertainties. 

Subject headings: methods: analytical — binaries: eclipsing — planets and satel- 
lites: general 



1. Introduction 

The transit of an exoplanet across the face of its parent star is an opportunity to learn 
a great deal about the planetary system. Photometric and spectroscopic observations reveal 
details about the planetary radius, mass, atmosphere, and orbit, as reviewed recently by 
Charbonneau et al. (2007). Transit light curves, in particular, bear information about the 
planetary and stellar radii, the orbital inclinati on, and the mean density of the star (Mandel 



& Agol 2000, ISeager & Mallen-Ornelasl (120031 ). Gimenez 2007). Additional planets in the 



system may be detected through gradual changes in the orbital parameters of the transiting 
planet (Miralda-Escude 2002, Heyl & Gladman 2007), or from a pattern of anomalies in a 
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collection of midtransit times (Holman & Murray 2005, Agol et al. 2005, Ford & Holman 
2007). 

In general, the parameters of a transiting system and their uncertainties must be esti- 
mated from the photometric data using numerical methods. For example, many investigators 
have used ^-minimization schemes such as AMOEBA or the Levenberg-Marquardt method, 
along with confidence levels determined by examining the appropriate surface of constant 
Ax 2 (see, e.g., Brown et al. 2001, Alonso et al. 2004) or by bootstrap methods (e.g., Sato et 
al. 2005, Winn et al. 2005). More recently it has become common to use Markov Chain Monte 
Carlo methods (e.g., Holman et al. 2006, Winn et al. 2007, Burke et al. 2007). However, 
even when numerical algorithms are required for precise answers, it is often useful to have 
analytic approximations for the parameters as well as their uncertainties and covariances. 

Analytic approximations can be useful for planning observations. For example, one may 
obtain quick answers to questions such as, for which systems can I expect to obtain the most 
precise measurement of the orbital inclination? Or, how many transit light curves will I 
need to gather with a particular telescope before the statistical error in the planetary radius 
is smaller than the systematic error? Now that nearly 50 transiting planets are known, we 
enjoy a situation in which a given night frequently offers more than one observable transit 
event. Analytic calculations can help one decide which target is more fruitfully observed, and 
are much simpler and quicker than the alternative of full numerical simulations. Analytic 
approximations are also useful for understanding the parameter degeneracies inherent in 
the model, and for constructing relatively uncorrelated parameter sets that will speed the 
convergence of optimization algorithms. Finally, analytic approximations are useful in order- 
of-magnitude estimates of the observability of subtle transit effects, such as transit timing 
variations, precession-induced changes in the transit duration, or the asymmetry in the 
ingress and egress durations due to a nonzero orbital eccentricity. 

Mandel & Agol (2005) and Gimenez (2007) have previously given analytic formulas for 
the received flux as a function of the relative separation of the planet and the star, but 
their aim was to provide highly accurate formulas, which are too complex for useful analytic 
estimates of uncertainties and covariances. Protopapas et al. (2007) provided an analytic 
and differentiable approximation to the transit light curve, but they were concerned with 
speeding up the process of searching for tran sits in large databases, rather than parameter 
estimation. ISeager fc Mallen-Ornelad (120031 ) presented an approximate model of a transit 
light curve with the desired level of simplicity, but did not provide analytic estimates of 
uncertainties and covariances. 

This paper is organized as follows. In §[2] we present a simple analytic model for a transit 
light curve, using a convenient and intuitive parameterization similar to that of Seager & 
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Mallen-Ornelas (2003). In § [31 we derive analytic approximations for the uncertainties and 
covariances of the basic parameters, and in § [4] we verify the accuracy of those approxima- 
tions through numerical tests. Our model assumes that the flux measurements are made 
continuously throughout the transit, and that stellar limb- darkening is negligible; in § 14.11 
and § 14.31 we check on the effects of relaxing these assumptions. In § [5] we derive some useful 
expressions for the uncertainties in some especially interesting or useful "derived" parame- 
ters, i.e., functions of the basic model parameters. In §[6] we present alternative parameter 
sets that are better suited to numerical algorithms for parameter estimation utilizing the 
analytic formalism given in § [31 We compare the correlations among parameters for various 
parameter sets that have been used in the transit literature. Finally, § 7 gives a summary 
of the key results. 



2. Linear approximation to the transit light curve 

Imagine a spherical star of radius with a uniform brightness and an unocculted flux 
/o- When a dark, opaque, spherical planet of radius R p is in front of the star, at a center-to- 
center sky-projected distance of zR±, the received stellar flux is F e (r, z, fo) = fo(l — X e (r, z)), 
where 

{0 1 + r < z 

\ (r 2 K + y ^ft±^g ) i_ r<z <i + r f (i) 
r 2 z < 1 — r 

with ki = cos _1 [(l - r 2 + z 2 )/2z] and k = cos -1 [(r 2 + z 2 - l)/2rz] (Mandel & Agol 2002). 
Geometrically, A e is the overlap area between two circles with radii 1 and r whose centers 
are z units apart. The approximation of uniform brightness (no limb darkening) is valid for 
mid-infrared bandpasses, which are increasingly being used for transit observations (see, e.g., 
Harrington et al. 2007, Knutson et al. 2007, Deming et al. 2007), and is a good approximation 
even for near-infrared and far-red bandpasses. We make this approximation throughout this 
paper, except in § 14.31 where we consider the effect of limb darkening. 

For a planet on a circular orbit, the relation between z and the time t is 

z(t) = aR~ l y/[sin n{t - t c )) 2 + [cos icos n{t-t c )) 2 (2) 

where a is the semimajor axis, % is the inclination angle, n = 27r/P is the angular frequency 
with period P, and t c is the transit midpoint (when z is smallest). 

The four "contact times" of the transit are the moments when the planetary disk and 
stellar disk are tangent. First contact (ti) occurs at the beginning of the transit, when the 
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disks are externally tangent. Second contact (tn) occurs next, when the disks are internally 
tangent. Third and fourth contacts (tni and try) are the moments of internal and external 
tangency, respectively, as the planetary disk leaves the stellar disk. The total transit duration 
is tiy — tj. The ingress phase is defined as the interval between t\ and tn, and likewise the 
egress phase is defined as the interval between tm and tjy. We also find it useful to define 
the ingress midpoint t ing = (ti + tu)/2 and the egress midpoint t egT = (tm + £rv)/2. 

Although Eqns. ([T]) and ([2]) give an exact solution, they are too complicated for an 
analytic error analysis. We make a few approximations to enable such an analysis. First, 
we assume the orbital period is large compared to transit duration, in which case Eqn. (jSJ) 
is well-approximated by 

z(t) = VW zr Q/W+^, (3) 

where, for a circular orbit, tq = R :k P/2ira = R+/na and b = acosi/R± is the normalized 
impact parameter. In this limit, the planet moves uniformly in a straight line across the 
stellar disk. Simple expressions may be derived for two characteristic timescales of the 
transit: 

W - ting = r (V(l + rf - V + ^ (I - rf - b^ = 2t VT^¥ + 0(r 2 ) (4) 
tu-ti = r (V (1 + r) 2 - 6 2 - y/(l - r) 2 - 6 2 ) = ^ -^== + 0(r 3 ). (5) 

It is easy to enlarge the discussion to include eccentric orbits, by replacing a by the planet- 
star distance at midtransit, and n by the angular frequency at midtransit: 

a(l-e 2 ) 

a — > , 

1 + e sin to 

n(l + esinu) 2 



n 



'1 -e 2 )i 



where e is the eccentricity, and uj is the argument of pericenter. Here, too, we approximate 
the planet's actual motion by uniform motion across the stellar disk, with a velocity equal 
to the actual veloci ty at midtransit. Methods for computing these quantities at midtran- 



sit are discussed by iMurray fc Dermottl (120001 ) , as we ll as recent t r ansit- specific studies by 



Barnes (2007), Burke (2008), Ford et al. (2008), and iGillon et all (120071 ). We redefine the 
parameters tq and b in this expanded scope as 

a cos i ( 1 — e 2 

(6) 



R± V 1 + e sin u 



„ 3 * (V^ZE. , . (7) 



an V 1 + e sin uj 



We do not restrict our discussion to circular orbits (e = 0) unless otherwise stated. 



Next, we replace the actual light curve with a model that is piecewise-linear in time, as 
illustrated in Figure [TJ Specifically, we define the parameters 



5 = f r 2 = f (R p /R* 



T 



r 

2r . 



(9) 
(10) 



and then we define our model light curve as 



F\t) 



fo-S 



fo-S \t- t c \ < T/2 - r/2 

*(|*-* c | -T/2 + T/2) T/2 - r/2 < \t - t c \ < T/2 + r/2 (11) 



fo 



\t-t c \ > T/2 + r/2 



We use the symbol F l to distinguish this piecewise-linear model (I for linear) from the 
exact uniform- source expression F e given by Eqns. and (T5]). The deviations between F l 
and F e occur near and during the ingress and egress phases. The approximation is most 
accurate in the limit of small r and b and is least accurate for grazing transits. As shown in 
Eqn. (JHJ), when r is small, r ~ tu — 1\ (the ingress or egress duration) and T ps t cgr — t ing (the 
total transit duration). Neithe r this piecewise-linear model nor the choice of parameters is 



new. 



Seager fe Mallen-Ornelasl (120031 ) also u sed a piecewi s e-line ar m odel, with dif f erent linear 



combinations of these parameters, and both lBurke et al.l (120071 ) and lBakos et al.l (120071 ) have 
employed parameterizations that are closely related to the parameters given above. What is 
specifically new to this paper is an analytic error and covariance analysis of this linear model, 
along with useful analytic expressions for errors in the physical parameters of the system. 
The "inverse" mapping from our parameterization to a more physical parameterization is 



a cos 2 



(Rp/R*) 2 

2 \ 2 



S/fo 



1 - e 



(12) 
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Fig. 1. — Comparison of the exact and piecewise-linear transit models, for the parameter 
choice r = 0.2, b = 0.5. The dashed line shows the exact uniform-source model F e , given by 
Eqn. (Q. The solid line shows the linear model F l , given by Eqn. (fTTl) . 
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3. Fisher information analysis 



Given a model F(t; {pi}) with independent variable t and a set of parameters {pi}, 
it is possible to estimate the covariance between parameters, Cov(pi,pj), that would be 
obtained by measuring F(t) with some specified cadence and precision. (Gould 2003 gives 
a pedagogical introduction to this technique.) Suppose we have iV data points taken at 
times tk spanning the entire transit event. The error in each data point is assumed to be a 
Gaussian random variable, with zero mean and standard deviation a^. Then the covariance 
between parameters {p{\ is 



where B is the zero-mean Gaussian-noise Fisher information matrix, which is calculated as 



Here, Bu is the inverse covariance matrix of the flux measurements. We assume the mea- 
surement errors are uncorrelated (i.e., we neglect "red noise"), in which case Bki — $ki&k ■ 
We further assume that the measurement errors are uniform in time with cr^ = a, giving 
Bki = <W~ 2 - 

In Table (JTJ, we compute the needed partial derivative^ of the piecewise-linear light 
curve F l , which has five parameters {pi} = {t c , r, T, 5, fo}. 



1 In computing these derivatives we have ignored the dependence of the piecewise boundaries in Table. (JTJ 
on the parameter values. The derivatives associated with those boundary changes are finite, and have 
a domain of measure zero in the limit of continuous sampling. Thus they do not affect our covariance 
calculation. 




(15) 




(16) 
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Totality Ingress/Egress Out of Transit 


£-F\t-{p m }) 
f T F% { Pm }) 
£F l (t; {p m }) 
j- s F\t-{p m }) 
£F l (t; {p m }) 




o -^{\t-t c \-l) 

-± 

- 1 Ml*-*d-f)-5 

1 1 1 



Table 1: Table of partial derivatives of the piecewise-linear light curve F, in the five parame- 
ters {p^ = {t c ,r,T,8,fo}. The intervals \t-t c \ < T/2-r/2, T/2-r/2 < \t-t c \ < T/2 + r/2, 
and \t — t c \ > T/2 + r/2 correspond to totality, ingress/egress, and out of transit respectively. 
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Fig. ([2]) shows the time dependence of the parameter derivatives, for a particular case. 
The time dependence of the parameter derivatives for the exact uniform-source model F e is 
also shown, for comparison, as are the numerical derivatives for limb-darkenened light curves. 
This comparison shows that the linear model captures the essential features of more realistic 
models, and in particular the symmetries. The most obvious problem with the linear model 
is that it gives a poor description of the r-derivative and the ^-derivative for the case of 
appreciable limb darkening, as discussed further in § 14.31 From Fig. <^ and Table flTJ we see 
that for the parameters T, r, and 5, the derivatives are symmetric about t = t c , while the 
derivative for the parameter t c is antisymmetric about t c . This implies that t c is uncorrelated 
with the other parameters. (This is also the case for the exact model, with or without limb 
darkening.) 
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Fig. 2. — Parameter derivatives, as a function of time, for the piecewise-linear model light 
curve F l (top row), the exact light curve for the case of zero limb darkening F e (second row), 
and for numerical limb-darkened light curves with a linear limb-darkening coefficient u = 0.2 
(third row) and u = 0.5 (bottom row). See § 14.31 for the definition of u. Typical scales are 
shown in the first row and are consistent in the following rows. 
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We suppose that the data points are sampled uniformly in time at a rate T = N/T to t, 
where the observations range from t = to to t = t Q + T tot and encompass the entire transit 
event. In the limit of large IV we may approximate the sums of Eqn. (fTSj) with time integrals, 
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to+Ttc 



to 



Using the derivatives from Table 
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dt. 



(17) 





-§ o 

1 -f 

rp T rp 

3 

—T Ttot / 



In what follows, it is useful to define some dimensionless variables: 

Q = 



TT-, 
a 



9 = r/T, 

V = T/{T tot -T-r). 



(19) 



The first of these variables, Q, is equal to the total signal-to-noise ratio of the transit in the 
limit r — > 0. The second variable, 8, is approximately the ratio of ingress (or egress) duration 
to the total transit duration. The third variable, 77, is approximately the ratio of the number 
of data points obtained during the transit to the number of data points obtained before or 
after the transit. Oftentimes, r and 9 are much smaller than unity, which will later enable 
us to derive simple expressions for the variances and covariances, but for the moment we 
consider the general case. 

Inverting B, we find the covariance matrix for the piecewise-linear model, 



Cov({t c , r, T, 5, f } , {t c , t, T, 5, f }) 


b - is] ^ 
b + A] e5T 

?l65T 



1 

~Q 2 



( -T 2 






\ 

[ry-^j^T 2 [rj + ^dST V 65T 

rj65T v 5 2 r]5 2 ) 



(20) 



The elements along the diagonal of the covariance matrix are variances, or squares of standard 
errors, a Pi = a/Cov^,^). 
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This result can be simplified for the case when many out-of-transit observations are 
obtained and rj — > 0. In this limit, Jq is known with negligible error, and we may assume 
fo = 1 without loss of generality In this case, 5 is the fractional transit depth, and the 
covariance matrix becomes 



/ e T 2 



Cov({t c , r, T, 5} , {t c , t, T, 5}) = ± 














61(6-50) T 2 
1-9 J 


e 2 T 2 
i-e 1 





e 2 rp2 

i-e 1 


e(2-e) T 2 
i-e 1 




i-e u± 


i-e u± 





i-e u± 

°—§T 

i-e u± 

i-e u 



\ 



■ (21) 



/ 



from which it is obvious that 9 is the key controlling parameter that deserves special atten- 
tion. Using Eqns. (Q and (jTOj) we may write 



9 



b 2 ' 



(22) 



Unless the transit is grazing, we have b < 1 — r, and 9 is restricted to the range [r , ^-1 . 
Fig. ([3]) shows the dependence of 9 on the impact parameter, for various choices of the transit 
depth. It is important to keep in mind that for b < 0.5, 9 is nearly equal to r and depends 
weakly on b. This implies that 9 is expected to be quite small for most transiting systems. 
For planetary orbits that are randomly oriented in space, the expected distribution of b is 
uniform, and hence we expect 9 < 0.3 for 90% of a random sample of transiting planets 
with R p < -RjuJE For this reason, in the following figures we use a logarithmic scale for 
9, to emphasize the small values. Fig. (j3J) shows the (suitably normalized) elements of the 
covariance matrix as a function of 9. 



2 In fact, the fraction of discovered systems with 9 < 0.3 may be even larger than 90%, because selection 
effects make it harder to detect grazing transits. 




Fig. 3. — Dependence of 9 = ^ on depth 5 = r 2 and normalized impact parameter b, for the 
cases r = 0.05 (solid line), r = 0.1 (dashed line), and r = 0.15 (dotted line). 
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Fig. 4. — Standard errors and covariances, as a function of 9 = t/T, for different choices of 
rj. The analytic expressions are given in Eqn. ( 1201) . The definitions of rj, 6, and Q are given 
in Eqn. (|T9l) . Solid line — = 0; Dashed line — i] = 0.5; Dotted line — r) — 1. 
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In the limits rj — > (errorless knowledge of /o) and 9 —> r (small impact parameter) 
the expressions for the standard errors are especially simple: 



Q^Ty/W, 



(23) 



In this regime, we have a clear hierarchy in the precision with which the time parameters 
are known, with a tc < < u T . 

To further quantify the degree of correlation among the parameters, we compute the 
correlation matrix, 



Corr({t c , r, T, S, /„} , {t c , r, T, 5, /„}) 

(\ 
1 



Cov(i,j) 



y/Cov(i,i)Cov(j,j) 


09-1)9 



\ 



v /(6-e(5-/3))(2-e(l-/3)) 



03-1)0 



03+1)0 
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(g-jVg 



6-9(5-/3) 



v / (6-e(5-/3) )(2-g(l-/3)) 
09+1)0 



og-iyg 



v /(/3+l)(2-e(l-/3)) V 2-0(1-/3) 



.(24) 



6-0(5-/3) 



6-0(5-/3) 



v /(/3+i)(2-0(i^)y 



1 
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/3+1 



/3+1 



2-0(1-/3) 

where we have defined (3 = 7/(1 — 0) to simplify the resulting expression. For 9 —>■ 0, all 
correlations with f vanish except for the correlation with 5. Due to the fact the correlation 
between 5 and /o is oc P 1 ^ 2 , it remains large even for fairly small (3. In the limit of r\ — ► 
(/? — > 0), we remove all correlations with /o and have the remaining correlations depending 
only on the ratio 9: 

( 1 
1 



lim Corr(- 

77->-0 











v /(6-50)(2-0) 



6-50 










V o o 



v /(6-50)(2-0) 



-50 



1 








1 / 



(25) 



Correlations with / decline with r\ as y/rj. 



In Fig. we have plotted the nonzero correlations as a function of 9 for a few choices 
of r). The special case of 77 — ^ is plotted in Fig. ([6]). In the 7/ — > limit, all correlations 
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are small (< 0.3) over a large region of the parameter space. Thus, our choice of parameters 
provides a weakly correlated set for all but grazing transits (6 ~ 1/2), as noted during the 
numerical analysis of particular systems by Burke et al. (2007) and Bakos et al. (2007). One 
naturally wonders whether a different choice of parameters would give even smaller (or even 
zero) correlations. In § [6] we present parameter sets that are essentially uncorrelated and 
have other desirable properties for numerical parameter estimation algorithms. 

The analytic formalism given in this section and more specifically the simple analytic 
covariance matrices in Eqns. (12UI I2T|) provide a toolbox with which to evaluate the statistical 
merits of any parameter set that can be written in terms of our parameters. In § [5] this 
technique is defined and applied to produce analytic formulas for variances, covariances and 
uncertainties in several interesting parameters. 




Fig. 5. — Correlations of the piecewise-linear model parameters, as a function of 9 = r/T 
for different choices of rj. Solid line — rj = 0; Dashed line — rj = 0.5; Dotted line — rj = 1. 
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Fig. 6. — Correlations of the piecewise-linear model parameters, as a function of 9 = r/T, 
for the case 77 — > (errorless knowledge of the out-of-transit flux). Solid line — Corr(r, T). 
Dashed line — Corr(r, S). Dotted line — Corr(T, S). 
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4. Accuracy of the covariance expressions 



Before investigating other parameter sets, it is necessary to examine the validity of 
Eqns. f l20| [2TI l24"t [251) when compared to similar quantities derived from more realistic 
transit light curve models. The utility of the covariance matrix in Eqn. (1201) depends on 
the accuracy of the integral approximation of Eqn. (|17|) . and on the fidelity with which the 
parameter dependences of the piecewise-linear model mimic the dependences of the exact 
uniform- source model. In this section we investigate these two issues. 



The case of a finite observing cadence, rather than continuous sampling, can be analyzed 
by evaluating the exact sums of Eqn. ([TBI) . Generally, given a sampling rate T, we expect 
the integral approximation in Eqn. (fTTj) to be valid to order (IY) -1 . In the 77 — ^ limit 
we may evaluate the exact sums, under the assumption of a uniform sampling rate, with 
data points occurring exactly at the start and end of the ingress (and egress) phases as well 
as at some intermediate times. This directly summed covariance, Cov sum , is related to the 
integral-approximation covariance Eqn. (1211) as 



where e = (IV) 1 . 

The quantity IV is approximately the number of data points obtained during ingress or 
egress. It is evident from Eqn. (I26p that for this sampling scheme only the variances of T and 
t along with their covariance are corrected. The corrections to the variances and covariance 
are 0(e 2 ) and 0(e) respectively. 

4.2. Comparison with covariances of the exact uniform-source model 

We tested the accuracy of the covariance matrix based on the piecewise-linear model by 
(1) performing a numerical Fisher analysis of the exact uniform-source model, and also (2) 
applying a Markov Chain Monte Carlo (MCMC) analysis of simulated data based on the 
exact uniform- source model. In both analyses, orbits are assumed to be circular. For the first 
task, we evaluated the analytic parameter derivatives of Eqn. ([1]), which are too cumbersome 



4.1. 
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(26) 
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to be worth reproducing here, and numerically integrated Eqn. (fTTj) to generate covariance 
matrices over a wide range of parameter choices. Fig. (j2J), in § 2, shows the parameter 
derivatives for the exact model, as well as the piecewise-linear model and some limb-darkened 
light curves. For the second task, idealized data was generated by adding Gaussian noise 
with standard deviation o / fo = 5 x 10~ 4 to Eqn. (pQ) sampled at V — 100 (in units of the 
characteristic timescale tq, Eqn. ([7])). With this sampling rate, approximately 50 samples 
occur during the ingress and egress phases. Approximately 10 4 links per parameter were 
generated with a Gibbs sampler and a Metropolis-Hasting jump-acceptance criterion. The 
jump-success fraction (the fraction of jumps in parameter space that are actually executed) 
was approximately 25% for all parameters. The effective length, defined as the ratio of 
the number of links to the correlation length (see the end of § [6] for the exact definition), 
was roughly 1000 — 2000 for the piecewise-linear mode l par ameter set. M ore details on the 



MCMC algorithm are given by iTegmark et al.l (120041 ) and iFordl (120051 ) . Standard errors 



were determined by computing the standard deviation of the resulting distribution for each 
parameter. The Fisher-information analysis should mirror the MCMC results, as long as the 
log-likelihood function is well approximated as quadratic near the mean (Gould 2003). 

The numerical Fisher analysis was performed for rj = and 0.05 < 9 < 1/2. In practice 
this was done by choosing r = 0.05 and varying b across the full range of impact parameters. 
(The numerical analysis confirmed that the suitably-normalized covariances vary only as a 
function of 9 = t/T, with the exception of slight (^-dependent positive offset in as that goes 
to zero as 5 goes to zero.) The MCMC analysis for rj = was accomplished by fixing the 
out of transit flux, f = 1, and varying the remaining parameters We chose r = 0.1 for the 
MCMC analysis. Fig. ([7]) shows all of the nonzero numerical correlation matrix elements, 
as a function of 9. The MCMC results, plotted as solid symbols, closely follow the curves 
resulting from the numerical Fisher analysis. Fig. ([8]) shows the nonzero numerical covariance 
matrix elements, also for the case 77 = 0. 

The correlations of the piecewise-linear model match the correlations of the exact model 
reasonably well, with the most significant deviations occurring only in the grazing limit, 
9 ~ 1/2. We have also confirmed that a similar level of agreement is obtained for nonzero rj, 
although for brevity those results are not shown here. We concluded from these tests that 
the errors in the analytic estimates of the uncertainties are generally small enough for the 
analytic error estimates derived from the piecewise-linear model to be useful. 
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Fig. 7. — Comparison of the non-zero correlation matrix elements for the exact light-curve 
model and the piecewise-linear model, as a function of 9 = r/T, for rj — > 0. Black curves: cor- 
relations for the piecewise-linear model. Gray curves: correlations for the exact uniform- 
source model. Black dots: correlations based on an MCMC analysis of simulated data with 
Gaussian noise (r = 0.1). 



-22 - 




Fig. 8. — Comparison of the covariance matrix elements for the exact uniform-source model, 
linear limb-darkened model, and the piecewise- linear model, as a function of 9 = t/T, for 
?7 — ^ 0. Black curves: covariances for the piecewise-linear model. Gray curves: covariances for 
the exact model with linear limb-darkening coefficient u = (solid) and u = 0.5 (dashed). 
Black dots: covariances as determined by a MCMC analysis of simulated data with Gaussian 
noise (u — and r = 0.1). The dimensionless number Q = \/TT5/a (see Eqn. [T9|) is 
approximately the signal-to-noise ratio of the transit. 
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4.3. The effects of limb darkening 



The piecewise-linear function of Eqn. ffTTj) was constructed as a model of a transit across 
a stellar disk of uniform brightness, with applications to far-red and infrared photometry in 
mind. At shorter wavelengths, the limb darkening of the star is important. How useful are 
the previously d erived results for this c ase, if at all? We used the limb-darkened light-curve 



models given by iMandel fc Agoll (120021 ) to answer this question. 



To simplify the analysis we adopted a "linear" limb-darkening law, in which the surface 
brightness profile of the star is 



I{z) 



l-u(l- VI - z 2 ^j (27) 



where u is the linear limb-darkening parameter. Claret (2000) finds values of u ranging from 
0.5-1.2 in UBVR for a range of main-sequence stars. Longer wavelength bands correspond 
to a smaller u for the same surface gravity and effective temperature. Solar values are u ~ 0.5 
in the Johnson R band and 0.2 in the K band. Fig. ([2]) of § 2 shows the time-dependence of 
the parameter derivatives of a linear limb-darkened light curve, for the two cases u = 0.2 and 
u = 0.5, to allow for comparison with the corresponding dependences of the piecewise-linear 
model and the exact model with no limb darkening. 

From the differences apparent in this plot, one would expect increased correlations 
(larger than our analytic formulas would predict) between the transit depth and the two 
timescales r and T. This is borne out by our numerical calculations of the covariance matrix 
elements, which are plotted in Figs. ( J51I91) . The analytic formulas underpredict the variances 
in 5 and r by a factor of a few, and they also severely underpredict the correlation between 
those parameters. 
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Fig. 9. — Comparison of the analytic correlations (black lines; Eqn. 12^1) numerically- 
calculated correlation matrix elements for a linear limb-darkened light curve (gray lines), 
as a function of 9 = t/T, for ?7 — ^ 0. Linestyles follow the conventions of Fig. (171) . 
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It is possible to improve the agreement with the analytic formulas by associating S 
with the minimum of the transit light curve, rather than the square of the radius ratio. 
Specifically, one replaces the definition of Eqn. (|12|) with the new definition 



„ 9 - 8 (VT^¥ - 1) u , s 

5 = f r 2 K - '—. 28 

J 9 -8m v ' 

For the previously- derived formulas to be valid, we must adopt a value for u based on 
other information about the parent star (its spectral energy distribution and spectral lines, 
luminosity, etc.) rather than determining u from the photometric data. Fig. ([TO]) shows the 
correlations resulting from this new association, for the case u = 0.5. Fig. ( JTTj) shows the 
improvement with this new association for the variance in S and the covariance between 6 
and r, for the case u = 0.5. While this new association improves on the agreement with 
the analytic covariances (particularly at low normalized impact parameter), a disadvantage 
is that we no longer have a closed-form mapping from {5, T, r} back to the more physical 
parameters {r, b, r }. 

It should be noted that there is evidence that linear limb darkenin g may not adequatel y 



fit high-qua li ty tra nsit light curves relative to higher order models (IBrown et al.l (120011 ) 



Southworthl (120081 )). A more complete analysis with arbitrary source surface brightness 



wou l d min imally include quadratic limb darkening but is outside the scope of this discussion. 



Pall (120081 ) completes a complementary analysis to this one of uncertainties in the quadratic 



limb darkening parameters themselves. 
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Fig. 10. — Comparison of correlation matrix elements for the piecewise-linear model (black 
curve) and a linear limb-darkened light curve (u = 0.5; gray curve), as a function of 9 = r/T. 
Here, the 5 parameter has been redefined as the minimum of the limb-darkened light curve, 
as approximated by Eqn. (1251) . Linestyles follow the conventions of Fig. (j7J). 
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Fig. 11. — Comparison of select covariance matrix elements for the piecewise-linear model 
(black curve) and a linear limb-darkened light curve (u = 0.5; gray curves), as a function of 
9 = r/T. The 5 parameter has been redefined as the minimum of the limb-darkened light 
curve, as approximated by Eqn. (128]) . in the solid gray curve. The dashed gray curve uses 
the initial S association, as defined in Eqn. ([8]). Linestyles follow the conventions of Fig. (JSJ). 
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5. Errors in derived quantities of interest in the absence of limb darkening 

The parameters {t c , r, T, 6, fo} are preferred mainly because they lead to simple analytic 
formulas for their uncertainties and covariances. The values of these parameters are also 
occasionally of direct interest. In particular, when planning observations, it is useful to 
know the transit duration, depth, and the predicted midtransit time. Of more direct scientific 
interest are the values of the "physical" parameters, such as the planetary and stellar radii, 
the orbital inclination, and the mean density of the star. Those latter parameters also offer 
clearer a priori expectations, such as a uniform distribution in cosi. 

For affine parameter transformations p i— > p' , we may transform the covariance matrix 
C via the Jacobian J = IP- as 

dp 

C = J T CJ . (29) 



Using Eqns. p2Hl4|) . we may calculate the Jacobian 



d{t c ,b 2 ,T$,r,f } 
d{t c ,r,T,6J } 
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(30) 



between the parameters of the piecewise-linear model and the more physical parameter set 
when limb darkening is negligible. Using this Jacobian, the transformed covariance matrix 
is 
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where we have ignored the unmodified covariance elements involving t c , and have kept only 
the leading-order terms in r in the //-dependent matrix. 
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The standard errors for other functions of the parameters, f{{pi}), can be found via 
error propagation, just as in Eqn. (129]) . 



The results for several interesting and useful functions, such as the mean densities of the 
star and planet, are given in Table ([2]). For brevity, the results are given in terms of the 
matrix elements of Eqn. (I3TI) . Simplified expressions for covariance matrix elements in the 
limit of 1] — > 0, 9 small (plentiful out-of-transit data) and negligible limb darkening are given 
in Table ©. 




(32) 



- 30 - 



f 1 1 1 Q nT 1 il T 

LlclIILlL^y 


V cLl IdllCt! I O LdllLlcLl LI IUi Lldl fcJU, 1 


1>I Ulco 


R p = rR+ 


R 2 p [Var(r)/r 2 + (log MjM ) 2 Var(x)] 


1 


P*/a = (7i/7 2 )27rro/P 


3 OH*/ ) Var ( r o)/ r o 




P P /a = ( 7 i/7 2 )27rr r/P 


(i? p /a) 2 [iVar(r 2 )/r 4 + Var(r)/r 2 ] 




= (72/7i)|acosi/-R*| 


±Var(& 2 )/6 2 




cos i 

= ( 7l 2 /7 2 3 )27rr |6| /P 


i cos 2 i [Var(6 2 )/6 4 + Cov(r 2 , 6 2 )/r 2 6 2 + Var(r 2 )/r 4 ] 




P* 

= (7 2 /7i) 3 (3/8G7r 2 )P/r 3 


fp 2 Var(r 2 )/r 4 




Pp 

= l2 (K^/r 3 sim) (P/2nGM^ 3 


p 2 [J Var(r 2 )/r 4 + 9Var(r)/r 2 + |Oov(r, r 2 )/rr 2 

il/ • / 7 \ 4t r / 7 2 \ 3 / ■ n \ 2n 111 2 \ / 2 

+ |(cosz/o) Var(o ) - |(cosz/o) Cov(o ,To)/tq 
-|(cosi/6) 2 Cov(6 2 , r)/r + Var(fQ/JY 2 ] 


2 


= ( 72 /7i) 3 P^/(2vrr 3 ) 


gl [|Var(r 2 )/r 4 + (log M*/M ) 2 Var(x)] 


1 


9 P 

= ( 1 l/ 1 2 )K*P/(2ixr 2 T 2 s\ni) 


(7 2 [Var(r 2 )/r 4 + 4Var(r)/r 2 + 2Cov(r, r 2 )/rr 2 
+ i(cosV6) 4 Var(6 2 ) - |(cos t/b) 2 Cov{b 2 , r 2 )/r 2 
-(cosz/&) 2 Cov(6 2 , r)/r + Var(iYj/iT 2 ] 


2 



Table 2: Table of transit quantities and associated variances, in terms of the matrix elements 
given in Eqn. ( |3TI) . We have assumed that both the orbital period, P, and stellar mass, 
M*, are known exactly. We have defined the noncircular-orbit parameters 71 = 1 + esina; 
and 72 = y/l — e 2 where e is the eccentricity and uo is the argument of pericenter (see § [2] 
for a discussion of eccentric orbits). Notes: (1) A mass-radius relation R± oc (M±/M & ) x 
is assumed; (2) We have assumed i > 80° in simplifying the inclination dependence in the 
variance. Quantities in bold are not determined by the transit model and must be provided 
from additional information. is the semi-amplitude of the source radial velocity. Terms 
have been arranged in order of relative importance with the largest in absolute magnitude 
coming first. Refer to Table ([3]) for matrix elements of Eqn. fl3Tl) for the case in which the 
planet is small, the out-of-transit flux is known precisely and limb darkening is negligible. 
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Table 3: Covariance matrix elements from Eqn. ( 1311) in the limit 77 — and 9 small for use 
in Table (J2]). These approximations are valid in the case in which the planet is small, the 
out-of-transit flux is known precisely and limb darkening is negligible. 
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6. Optimizing parameter sets for fitting data with small limb darkening 

The parameter set {t c , r, T, S, fo} has the virtues of simplicity and weak correlation 
over most of the physical parameter space. However, when performing numerical analyses of 
actual data, the virtue of simplicity may not be as important as the virtue of low correlation, 
which usually leads to faster and more robust convergence. To take one example, lower 
correlations among the parameters result in reduced correlation lengths for Monte Carlo 
Markov Chains, and faster convergence to the desired a posteriori pro bability distributions , 



and can obviate the need for numerical Principal Component Analysis (ITegmark et al.ll2004j ) . 
In Fig. (TT2"j) . we compare the degree of correlations for various parameter sets that have been 
used in the literature on transit photometry. Of note is the high degree of correlations among 
the "physical" parameter set {R*/a, Rp/a, b}, which is a poor choice from the point of view 
of computational speed. 
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Fig. 12. — Comparison of correlations for various parameter sets that have been used in the 
literature. The correlations were derived from the piecewise-linear model (Eqn. |2"T|) assuming 
71 = 0. (a) Parame ters {b 2 ,TQ,r}. (b) {R±/a = nr ,R p /a = nr r,b 2 }. (c) {2/T, b 2 ,r} (e.g., 
Bakos et al.l (120071 )). (d) {T, r, 5}, the set introduced in this paper. 
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Nevertheless, one advantage of casting the model in terms of physical parameters is that 
the a priori expectations for those parameters are more easily expressed, such as a uniform 
distribution in b. The determinant of the Jacobian given by Eqn. (12"§1) . |J|, is also useful 
in translating a p rior i probability distributions from one parameter set to the other [see 



Burke et al.l (120071 ) or iFordl (120061 ) for an example of how this is done in practice]. For the 



case of the parameter set {t c , r, T, S, fo}, we may use the Jacobian, Eqn. (I30I) . to convert a 
priori probability distributions via 



p(t c , r, T, 5, f Q )dt c dr dT d5 df = p(t c , b , r , r, /, 



1 



-dt c db 2 c/tq dr dfo 
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where we have remeasured the phase space volume via the determinant, 
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(34) 



One may use this expression to enforce a uniform prior in b, for example, by weighting 
the likelihood function as shown in Eqn. (|33|) . However, there is a practical difficulty due 
to the singularity at b = 0. One way to understand the singularity is to note that uniform 
distributions in r, T lead to a nearly uniform distribution in 9 = t/T, which highly disfavors 
b = 0; in order to enforce a uniform distribution in b, the prior must diverge at low b. Fig. fl3]) 
graphically captures the steep variation for small b with 9. Consider, instead, the parameter 
set {t c ,b,T,r = ^/5/JoJo} where, from Eqn. (JSJ), b 2 = 1 - rT/r. We may calculate the 
determinant of the Jacobian (not reproduced here) 



d{t c ,b,T,rJ } 



d{t c ,T,T,5,f } 
Combining this result with Eqn. (133]) 
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p(t c ,b,T,r,fo)dt c dbdTdrdf = p(t c ,b,r ,r, f 

= p(t c ,b,T ,r,f 



1 T 

1 - b 2 4ro~ 
1 



dt c db dr dr df 



dt c db dro dr dfo- (36) 



The singularity at b = has been removed with this parameter choice. There is a singularity 
at b = 1 instead, which is only relevant for near-grazing transits, and is not as strong of a 
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singularity because of the square root. We confirm that this parameter set also enjoys weak 
correlations, as shown in Fig. ffl3]) . and therefore this set is a reasonable choice for numerical 
parameter-estimation algorithms. The merits of other parameter sets, from the standpoint 
of correlation and a priori likelihoods, may be weighed in a similar fashion, using the simple 
analytic covariance matrix of Eqn. (1201) . and the appropriate transformation Jacobian, in 
combination with Eqn. ( |29|) . 
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Fig. 13. — Correlations for the parameter set {b,T,r}. The correlations were derived from 
the piecewise-linear model (Eqn. [21]) assuming r\ = 0. 
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If the issues associated with the transformation of priors are ignored (i.e. if the data are 
of such quality that the results will depend negligibly on the priors), we can give essentially 
uncorrelated parameter sets. Consider, for example, the parameter set {t c , S e = 6/r,T, A = 
ST}. The new parameter S e is the magnitude of the slope of the light curve during the 
ingress and egress phases, and the new parameter A is the area of the trapezoid defined 
by the transit portion of the light curve (i.e., the time integral of the flux decrement). For 
simplicity we assume 77 = and fix fo = 1. The transformed correlation (Eqn. ([25]) ) is found 
via the transformation Jacobian, Eqn. (I2"9"j) as 



Corr({t c ,S e ,T,A} , {t c , S e , T, A}) 



( 1 
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The determinant of the transformation Jacobian (for use with Eqn. (133]) ) is given as 



d{t c ,S e ,T,A} 
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(38) 



With this new parameter set, the only nonzero correlation is between T and A, and this 
correlation is < 0.3 even for grazing transits (see Fig. [Til) . We have found that these param- 
eters provide a nearly optimal set for data fitting when little is known at the outset about 
the impact parameter of the transit. 

It is possible to do even better when the impact parameter is known at least roughly. 
Consider the parameter set {t c , S e , II = T5°, 5} where S e is the slope of ingress, and 6 is a 
constant (whose chosen value will be discussed momentarily). The new parameter II has no 
simple physical interpretation. We again assume rj = and /o = 1. The correlation matrix 
in this case is 



Corr({t c ,S e ,n,5} ,{t c ,S e ,U,S}) 
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The determinant of the transformation Jacobian is given as 



(39) 
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(40) 
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With this choice, the only nonzero correlation is between II and 5. If the constant 9 is 
chosen to be approximately equal to 9, then this sole correlation may be nullified. Thus, if 9 
is known even approximately at the outset of data fitting — from visual inspection of a light 
curve, or from the approximation 9 ~ r valid for small planets on non-grazing trajectories — 
a parameter set with essentially zero correlation is immediately available. As an example, 
Fig. (TH1) shows the correlation between IT and 5 as a function of 9, for the choice 9 = 0.1, 
which has a null at 9 = 0.1 as expected. 

The utility of this parameter set is not lost if 9 cannot be confidently specified when 
used with Markov chain Monte Carlo parameter estimation codes. At each chain step i, the 
next candidate state can be drawn from the candidate transition probability distribution 
function generated by the above parameter set with 9 = 9i-\. Thus, the Markov chain will 
explore the parameter space moving along principal axes at each chain step. Additionally, 
allowing the candidate transition function to vary as the Markov chain explores parameter 
space may prove useful for low S/N data sets. 
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Fig. 14. — Comparison of the correlations among the parameters, for the set {5, T, r} (black 
lines), the set {S e ,T, A = T5} (dashed-dot gray line) and the set {S e , II = T5 9 ,5} (solid 
gray line) for the case 9 = 0.1. For the latter set, the only nonzero correlation is between II 
and S e , which vanishes at 9 = 0.1. 
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As a concrete example of the effectiveness of uncorrelated parameters, we apply the 
MCMC algorithm to simulated data. For a given choice of the parameter set, we generate 
chains with a fixed jump-success fraction, and calculate the resulting autocorrelations of the 
Markov chain. For a particular parameter p (with value Pi at chain step i), the autocorrela- 
tion a at a given chain step j is defined as 

3 ~ (p 2 >-(p*> 2 { } 



where the averages refer to the averages over the whole chain ( Tegmark et al. 20041 ). The 



correlation length of the chain is the number of steps N that are required before the au- 
tocorrelation drops below 0.5. The total chain length divided by the correlation length is 
referred to as the effective length of a chain. The effective chain length is approximately the 
number of independent samples, which quantifies the degree of convergence of the algorithm. 
A lower correlation length, for the same total chain length, gives a more accurate final dis- 
tribution. This autocorrelation analysis was performed for both the "physical" parameter 
set {t c , b 2 , Tq, r 2 } as well as the parameter sets {t c , r, T, 5} and {t c , b, T, r}, with 77 = in all 
cases (i.e., plentiful out-of-transit data). The MCMC was executed as detailed in § 14.21 with 
a fixed jump rate ~50% for all parameter chains. (In practice this was achieved by adjusting 
the size of the Gaussian random perturbation that was added to each parameter at each 
trial step.) By choosing either the parameter set {t c ,r, T, 5} or {t c ,b, T, r}, the correlation 
lengths are reduced by a factor of approximately 150. By using the minimally-correlated 
parameter set {t c , S e , T, A}, the correlation lengths are reduced by an additional factor of 

To completely eliminate the correlations between parameters, one can diagonalize the 
symmetric covariance matrix, Eqn. (157)1 . a nd find the linear co mbinations of parameters that 



eliminates correlations. This was done by iBurke et al.l (120071 ) for the particular case of the 



transiting planet XO-2b. Analytic expressions for the eigenvectors are available because there 
are only two entangled parameters. However, these eigenvectors are linear combinations 
of local parameter values; they do not constitute a global transformation rendering the 
covariance diagonal. Thus, this procedure is useful for numerical analysis of a particular 
system, although not for analytic insights. 



7. Summary 

We have presented formulas for uncertainties and covariances for a collection of param- 
eters describing the light curve of an exoplanet transiting a star with uniform brightness. 
These covariances, given in Eqns. (T20| l3Tj) . are derived using a Fisher information analysis 
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of a linear representation of the transit light curve. The key inputs are the uncertainty in 
each measurement of the relative flux, and the sampling rate. We have verified the accuracy 
of the variance and covariance estimates derived from the piecewise-linear light curve with 
a numerical Fisher analysis of a more realistic (nonlinear) light-curve model, and with a 
Markov Chain Monte Carlo analysis of idealized data. 

We focused on a particular parameterization of this piecewise-linear light curve that we 
believe to be most useful. The parameters are the midtransit time (t c ), the out-of-transit 
flux (/o), the flux decrement during the full phase of the transit (S), the duration of ingress 
or egress (r), and the duration between the midpoint of ingress and the midpoint of egress 
(T). This set is observationally intuitive and gives simple analytic formulas for variances 
and covariances. The exact parameter definitions are provided in Eqns. ([HI EH HUj) in terms 
of normalized impact parameter, stellar and planetary radii, semi-major axis and orbital 
period. Inverse mappings to more physical parameters are provided in Eqns. (fl2| [T3| [T4|) . 
The analytic covariance matrix is given in Eqn. (|20|) and the analytic correlation matrix 
is given in Eqn. (1241) . Some quick-and-dirty (but still rather accurate) expressions for the 
parameter uncertainties, for the case in which the planet is small, the out-of-transit flux is 
known precisely and limb darkening is negligible, are given as 

a T « Q- l TV2§, 
as ~ Q- l 5 

where 9 = r/T is the ratio of the ingress or egress duration to the total duration, and 
Q = \/YT^ is the total signal-to-noise ratio of the transit in the small-planet limit (see 
Eqn. Mj. 

We investigated the applicability of these results to a limb darkened brightness profile, 
in which the true light curve is not as well-described by a piecewise-linear function. We 
found that the analytic formulas underestimate some of the variances and covariances by a 
factor of a few, for a typical degree of limb darkening at optical wavelengths. Significant 
improvements to covariance estimates in the limb darkened case may be made by redefining 
the depth parameter as a function of darkening coefficient and impact parameter as in 
Eqn. (128]) . Unfortunately, no closed-form mapping to more physical parameters exists with 
this choice, and therefore most of the appeal of the analytic treatment is lost. 

Quantities that are derived in part or in whole from the transit light curve (such as stellar 
mean density or exoplanet surface gravity) are provided in terms of the suggested parameter 
set. In Table fl2]), uncertainties propagated from the covariance estimates for these quantities 
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are provided with simple analytic formulas. In Table [3l covariance elements relevant to the 
uncertainties in Table [2] are given for the case in which the planet is small and the out-of- 
transit flux is known precisely. This allows the uncertainty in a given physical parameter 
to be predicted in advance of any data, bypassing the need for time-consuming simulations. 
For transit surveys, these formulas may also be useful in giving closed-form expressions for 
the expected distributions for some of the key properties of a sample of transiting planets. 

In §0 with the tools provided, we approach the question of what parameter sets are best 
suited to numerical parameter estimation codes. This question depends both on the level of 
parameter correlation and the behavior of any a priori likelihood functions. We advocated a 
parameter set that has the virtue of both weak correlation and essentially uniform a priori 
expectations: specifically, the parameters are the midtransit time, the out-of-transit flux, 
the ratio of planetary to stellar radii (Ftp/Ft*), the normalized impact parameter, and the 
duration between the midpoint of ingress and the midpoint of egress. Fig. (fT3l) graphically 
describes the parameter correlations while Eqn. (136]) gives the a priori probability distri- 
bution. Finally, two parameter choices are given that are less intuitive than the suggested 
set but that provide smaller correlations, depending on information that may be inferred 
or guessed prior to analysis. Correlations may be tuned to zero with the second parameter 
choice for a non-grazing transit and an estimate of R p / R*. The resulting correlation matrices 
for both parameter choices are given in Eqns. fl37l [39]) . Lower correlations relate directly to 
more efficient data fitting, as demonstrated by reduced correlation lengths with a Markov 
Chain Monte Carlo method. 
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